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Abstract. We investigate the effects of rotation on the behavior of the helium burning shell source in accreting carbon-oxygen 
white dwarfs, in the context of the single degenerate Chandrasekhar mass progenitor scenario for Type la supemovae (SNe la). 
We model the evolution of helium accreting white dwarfs of initially 1 M©, assuming four different constant accretion rates 
(2, 3, 5 and 10x10"^ M0/yr). In a one-dimensional approximation, we compute the mass accretion and subsequent nuclear 
fusion of helium into carbon and oxygen, as well as angular momentum accretion, angular momentum transport inside the 
white dwarf, and rotationally induced chemical mixing. Our models show two major effects of rotation: a) The helium burning 
nuclear shell source in the rotating models is much more stable than in corresponding non-rotating models — which increases 
the likelihood of accreting white dwarfs to reach the stage of central carbon ignition. This effect is mainly due to rotationally 
induced mixing at the CO/He interface which widens the shell source, and due to the centrifugal force lowering the density 
and degeneracy at the shell source location, b) The C/O-ratio in the layers which experience helium shell burning — which 
may affect the energy of a SN la explosion — is strongly decreased by the rotationally induced mixing of a-particles into the 
carbon-rich layers. We discuss implications of our results for the evolution of SNe la progenitors. 
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1. Introduction 

Accreting white dwarfs in close binary systems are consid- 
ered to be responsible for such important astrophysical phe- 
nomena as nova explosions, cataclysmic variables and super- 
soft X-ray sources. Especially, there is general agreement that 
accreting CO white dwarfs are progenitors of Type la super- 
novae (SNe la) (e.g., Woosley & Weaver ,1986. Branch et 
al. [19531 Whee ler fPy^ No moto et al. nWTl Hillebrandt & 
Niemever l2000l Livio l200ll l. 

The fairly homogeneous properties of light curves and 
spectra of SNe la have made them a valuable distance indicator 
to determine cosmological parameters (e.g, Hamuy et al. I996I 
Branch . 1998. Leibundgut .200 1 J . One of the most intriguing re- 
sults obtained from SNe la data may be the recent suggestion of 
a finite cosmological constant, which implies that the universe 
is currently accelerating its expansion (Perlmutter et al. .1999 
Riess et al.l2000b. 



The observational properties of SNe la demand a ther- 
monuclear explosion of a CO white dwarf as their origin. 
However, despite the far-reaching observational achievements, 
self-consistent progenitor evolution model calculations which 
explain the existence of SNe la are still lacking (but see Yoon 
& Langer l2003> . Models for the evolution of SN la progenitors 
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studied until now can be divided into two groups, one where 
the white dwarf explosion is triggered close to the center, and 
the other where the white dwarf explosion is triggered in the 
helium shell (see Livio l20()T1 for a review). 

Two scenarios exist for close-to-center ignition. One is the 
double degenerate scenario, in which two CO white dwarfs in 
a binary system merge due to the angular momentum loss by 
gravitational wave radiation. If the system mass exceeds the 
Chandrasekhar mass, carbon must ignite, with the possible out- 
come of a thermonuclear explosion. Although there is observa- 
tional evidence for double white dwarf systems (e.g., Koester 
et al. 1200 H Karl et al. 120031) . including one sdB star plus a 
CO white dwarf system which might result in a double white 
dwarf system having a total mass beyond the Chandrasekhar 
limit (Maxted et al. 2000 1, this scenario is criticized mainly by 
the fact that the merging of two white dwarfs is likely to end in 
an accretion-induced collapse rather than a SN explosion (Saio 
& NomotOil985,il998|l. 

The currently favored scenario for close-to-center ignition 
is the single degenerate scenario, in which a CO white dwarf in 
a close binary system accretes hydrogen or helium rich matter 
from a non-degenerate companion at a relatively high accretion 
rate {^ lO^^Mo/yr). The super-soft X-ray sources, which are 
believed to consist of a non-degenerate star and a white dwarf, 
are promising observed counterparts of this scenario. The idea 
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of steady nuclear burning due to a high accretion rate fits well to 
these sources (e.g., Li & van den Heuvei n"997l Kahabka & van 
den Heuvel 1997 Langer et al. 2000; Yoon & Langer 2003 1. 
Symbiotic systems (e.g., Hachisu et al..l999 1 and recurrent no- 
vae (e.g., Hachisu & Kato 20(5J Thoroughgood et al. 2001 1 
may also belong to this scenario. However, it is still not clear 
whether this scenario can explain the observed SN la frequency 
since the accretion rates which allow steady burning are lim- 
ited to a very narrow range. Furthermore, helium shell burning 
is generally found to be unstable: strong helium shell flashes 
might lead to strong mass loss by optically thick winds, or to a 
strong expansion of the white dwarf envelope to a red-giant 
phase which will result in a possible merging of the binary 
(e.g. Nomoto rT^82al Hachisu, Kato & Nomoto lD^ Cassisi 
et al. [T998l Kato & Hachisu[T999| Langer et al. 2002^. 

In the so called 'Edge Lit Detonation' models, corre- 
sponding to sub-Chandrasekhar mass models, helium ignites 
at the bottom of a helium layer of a few tenth of a so- 
lar mass, accumulated on a CO white dwarf at a relatively 
low accretion rate ( ^ IQ- ^Mp/yr) (e.g. No moto [T982bl 
Limongi & Tornambe 1991 Woosley & Weaver 119941 Livne 
& Arnett .l995j . These models are favored in terms of statis- 
tics (e.g. Iben & Tutukov [T99T1 Tout et al. |2001l but disfa- 
vored by the fact that the spectra and light curves produced by 
these models are not in agreement with observations (Hoflich et 
al. [T996l Nugent et al.ri997 Pinto et al.|200D- Yoon & Langer 



J200 4b) showed that a helium detonation in such models may 
be avoided altogether due to a strong shear layer inside an ac- 
creting white dwarf. 

Most previous studies of accreting white dwarfs as SN la 
progenitors have assumed spherical symmetry. However, in a 
close binary system which contains a white dwarf, an accre- 
tion disk is formed around the white dwarf. The matter coming 
from the accretion disk which is incorporated into the white 
dwarf carries high specific angular momentum, i.e. up to the 
value corresponding to Keplerian rotation at the white dwarf 
equator. Consequently, not only the accreted layer may be spin- 
ning very rapidly, but also the rest of the white dwarf, as angu- 
lar momentum may be transported into the stellar interior (e.g. 
Ritter .1985. Langer et al. .2000j . 

Recent polarization observations provide evidence for the 
aspherical nature of SNe la (SN 1999by, Howell et al.'200r 
SN 2001el, Wang et al 2003 1. A plausible explanation for this 
asphericity may be a rapidly rotating progenitor, while dumpi- 
ness of supernova ejecta or binarity of the progenitors could be 
another possible cause for the observed polarization (Kasen et 
al. |2003 1. The observation that white dwarfs in Cataclysmic 
Variables rotate much faster than isolated ones (Sion 1999 
Starrfield 2003 1 supports the idea that accreting white dwarfs 
are indeed spun up. In this picture, differential rotation in- 
side the accreting white dwarf is inevitable, unless the angular 
momentum transport in the white dwarf is extremely efficient 
such that rigid body rotation may be maintained throughout 
the white dwarf interior (Yoon & Langer "2004a'i. Rotationally 
induced instabilities such as the shear instability may lead to 
chemical mixing between the CO white dwarf material and 
the accreted matter. This effect has been discussed by sev- 
eral authors mainly with respect to the chemical enrichment 



of CNO and other heavier elements observed in nova ejecta 
(e.g., Kippenhahn & Thomas 1978 MacDonald [T983l Livio & 
Truran lTggTl Fuiimoto fT988l Fujimoto & Iben lT997t . 

In this paper, we investigate the effects of rotation on the 
evolution of accreting white dwarfs in relation to the single 
degenerate Chandrasekhar mass scenario for SNe la progeni- 
tors. For this purpose, we carry out numerical simulations of 
helium burning white dwarfs considering high accretion rates 
(> lO^^MQ/yr) which allow steady shell burning for hy- 
drogen or helium. Our discussions in this study are focused 
on the rotationally induced chemical mixing and their conse- 
quences for the behavior of the helium burning shell source. 
Comprehensive discussions on the rotational induced hydrody- 
namic instabilities and their role in the redistribution of angu- 
lar momentum in accreting white dwarfs are given in a sepa- 
rate paper (Yoon & Langer l2004a> . Effects of rotation on sub- 
Chandrasekhar mass models, i.e., on accreting white dwarfs 
with relatively low accretion rates (^ 10~^ Mo/yr), are dis- 
cussed in Yoon & Langer (■2004bJ . Our computational method 
is introduced in Sect. 2, and the initial model and physical 
assumptions are discussed in Sect. 3. In Sect. 4, we present 
the evolution of various physical quantities of our white dwarf 
models. We discuss rotationally induced chemical mixing and 
the stability of shell burning in Sect. 5 and 6 respectively. 
Finally, we discuss our results in Sect. 7. 



2. Method 

The numerical models have been computed with a hydrody- 
namic stellar evolution code (cf, Langer 1998> Heger, Langer 
& Woosley |2000j. Opacities are taken from Iglesias & Rogers 
(II996I I. Changes in the chemical composition are computed us- 
ing a nuclear network with more than 60 nuclear reactions, of 
which 3 1 reactions are relevant to helium burning such as 3a, 
i2C(a, 7)i^0(a, -f)^°Ne and 2"Ne(a, -ff^Mg{a, 7)28Sireac- 
tions. For the ^^C(a, 7)^^0 reaction rate, which has an impor- 



tant role in our study, we use the rate of Caughlan et al. (II985> 
reduced by a factor of 0.63, as suggested by Weaver & Woosley 
(19931. 

The effect of the centrifugal force on the stellar struc- 
ture, and rotationally induced transport of angular momen- 
tum and chemical species are treated in a one-dimensional 
approximation (Heger et al. 2000 1. We consider rotationally 
induced transport processes as Eddington-Sweet circulations, 
the Goldreich-Schubert-Fricke instability, and the dynamical 
and secular shear instability, by solving a non-linear diffusion 
equation, with diffusion coefficients for each of the instabili- 
ties computed as in Heger et al. ( 2000 1 and Yoon & Langer 
(I2004a> . The viscous dissipation of the rotational energy is also 
considered as described by Yoon & Langer (^2004aV More de- 
tails about the numerical methods are given in Yoon & Langer 
(l2004il . 

A comparison of our method with previously used ones 
for rotating white dwarf models is already given in Yoon & 
Langer (|2)04ai in detail. Here we give some additional com- 
ments on the rotation physics in relation to helium accreting 
white dwarfs. 
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Our approach to include rotation is based on the method 



given by Kippenhahn & Thomas J1970> . with the mass shells 
of a star defined by isobars instead of spherical shells. In de- 
scribing effects of rotationally induced chemical mixing, the 
chemical abundances are assumed to be constant on isobars. 
This is justified by theoretical achievements of Zahn ( 1975 1, 
Chaboyer & Zahn (|T992j and Zahn (1992 1, who showed 
that anisotropic turbulence acts much more efficiently in the 
horizontal than in the vertical direction. Shellular rotation 
and chemical homogeneity on isobars can be thus achieved 
(Meynet & Maeder 1997 ), enabling us to use a one dimensional 
approximation (see also Maeder l2003l . 

Our assumption of shellular rotation may be question- 
able for helium accreting white dwarfs, on two accounts. 
Kippenhahn & Thomas J1978> argued that the accreted mat- 
ter cannot spread over the whole surface of white dwarf but 
rather forms an accretion belt, since the angular momentum 
transport time scale based on the electron viscosity is longer 
than the accretion time scale. However, calculations by Sparks 
& Kutter c1987j failed in obtaining classical nova outbursts 
when the prescription of accretion belts was employed. Livio 
& Truran (1987 1 argued that if a turbulent viscosity is used, 
the angular momentum transport becomes much more efficient 
than discussed by Kippenhahn & Thomas. They concluded that 
the accreted matter can spread over the entire star with the 
consequence that mass accretion onto white dwarfs may re- 
sult in nova outbursts, which was not possible in the accre- 
tion belt prescription due to the too strong centrifugal force at 
the equator. A local stability analysis by MacDonald (1983) 
and a theoretical study by Fujimoto ii988t also concluded that 
the accreted matter will not remain confined to the equatorial 
belt due to efficient angular momentum transport toward the 
poles. We therefore conclude that the assumption of shellular 
rotation employed in our study may still be well-grounded, at 
least in the non-degenerate helium envelope, although multi- 
dimensional simulations are ultimately required for its justifi- 
cation (cf. Hujeirat 1995). 

Secondly, as degenerate matter is barotropic, perturbations 
which create dj /dz ^ 0, with j as specific angular mo- 
mentum and z being the distance to the equatorial plane, are 
eliminated on a short time scale. Thus, cylindrical rotation 
may be enforced in the inner CO core (e.g. Kippenhahn & 
Mollenhoff 1 1 9741 Durisen [l977> . This means that accreting 
white dwarfs may consist of two different rotation laws: shel- 
lular in the non-degenerate envelope and cylindrical rotation 
in the degenerate core. However, our numerical models can 
still represent the cylindrically rotating degenerate inner core to 
some degree, since most of the total angular momentum is con- 
fined at the equatorial plane in both the shellular and the cylin- 
drical cases. It is harder to estimate the angular momentum 
transport efficiency in the semidegenerate transition layer, as 
non-barotropic perturbations in the degenerate core may induce 
dynamical meridional flows (Kippenhahn & Mollenhoff" 1 974" 
Miiller 2003, private communication). Since it is poorly known 
how these meridional flows may affect the angular momentum 
redistribution, we have to leave the detailed investigation of this 
effect to future work. 



Finally, the accuracy of our results for very rapid rotation 
is limited due to the one dimensional approximation. In com- 
puting the effective gravitational potential in a rotating star, the 
potential is expanded in terms of spherical harmonics, of which 
we only consider terms up to the second order (Kippenhahn & 
Thomas [1 9701 Endal & Sofia 1976 1. This method can repro- 
duce the shapes of rigidly rotating polytropes accurately up to 
a rotation rate of about 60% of the critical value, correspond- 
ing to correction factors of fp ~ 0.75 and /t — 0.95 in the 
stellar structure equations (Fliegnerll993' Heger et al. I2000> . 
We therefore limit these factors to the quoted values, with the 
consequence that we underestimate the effect of the centrifu- 
gal force in layers which rotate more rapidly than about 60% 
critical. We shall see below that this affects only a very small 
fraction of our stellar models in most cases. However, as the 
surface layers are always very close to critical rotation, the stel- 
lar radius of our models may be somewhat underestimated. A 
better description of the equilibrium structure of rapidly rotat- 
ing stars which may deviate significantly from spheres requires 
multi-dimensional techniques such as the self-consistent field 
method (Ostriker & Mark 1 1 9681 Hachisu 1986 1. In addition, 
since our models are rotationally symmetric, three dimensional 
effects such as triaxial deformation can not be described with 
our numerical code. 



3. Physical assumptions 

We begin with a hot carbon-oxygen white dwarf of 
log is/L© — 3.379 in order to avoid the numerical difficulty of 
a strong initial helium shell flash. The initial mass of our model 
is 0.998 M0. The central temperature and density are initially 
such that Tc = 2.288 x 10** K and p^ = 1.64 x 10^ g cm-^, 
respectively. Although in reality white dwarfs may be much 
colder when mass accretion starts, the thermal evolution of the 
white dwarf core is not influenced by the different initial con- 
ditions due to the self heating of white dwarfs (cf. Yoon & 
Langer2 003i . Since observations (Heber et al. lI997l Koester et 
al. 1998 Kawaler'2003^ and stellar evolution models (Langer 
et al. 1999 1 find small spin rates of isolated white dwarfs, we 
adopt a slow and rigid rotation in our initial model: the initial 
rotation velocity at the white dwarf equator is set to 1 .0 km s~ ^ . 

We consider 4 different accretion rates: 2 x 10^^, 3 x 10^^, 
5 X 10^^ and 10^^ Mo/yr. The accreted matter is assumed to 
be helium-rich with Y ~ 0.98 and Z — 0.02. The abundances 
for heavier elements are chosen to be the same as in a helium 
core which results from core hydrogen burning in a main se- 
quence star through the CNO cycle. 

The specific angular momentum of the accreted matter 
is assumed to have the local Keplerian value, as most white 
dwarfs with a non-degenerate companion are believed to ac- 
crete matter through an accretion disk. However, analytic con- 
siderations can show this to cause an angular momentum 
problem. Assuming rigid rotation and a usual white dwarf 
mass-radius relation, the white dwarf reaches overcritical ro- 
tation much before reaching the Chandrasekhar limit (Livio & 
Pringle"1998" Langer et al. "2000*1. For finite angular momentum 
redistribution times, as realized by our preliminary calculations 
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Table 1. Properties of the computed models. The first column denotes the sequence number, where 'NR', 'R' and 'RT' denote the 
non-rotating models and rotating models with /c = 1/30 and 1/3 respectively. The second column gives the accretion rate in the 
unit of 10^^ M0/yr. The third column denotes the factor f^ for the chemical mixing efficiency in rotating models. Columns 4, 
5, 6 and 7 give mass, central density, central temperature and -Erot/|£^giav| of the last computed model of each evolutionary 
sequence. The eighth column gives the amount of accreted material until the onset of the helium shell instability. The column 
with LEdd indicates whether the last computed model reached the Eddington limit or not. The last column gives the total number 
of the computed models for each sequence. 



No. 



M 
[10-^ Me/yr] 



/c 



[M0] 



Pcjlast 

[10** g/cm^] 



-*- c,last 

[10* K] 



(-Brot/I-Bgrav Dlast 



AM 
[Mq] 



LEdd No. of Models 



NRl 


2.0 


- 


1.0026 


0.25 


1.45 


- 


0.00 


Yes 


4342 


NR2 


3.0 


- 


1.0042 


0.25 


1.45 


- 


0.00 


Yes 


5164 


NR3 


5.0 


- 


1.0130 


0.28 


1.37 


- 


0.00 


Yes 


10324 


NR4 


10 


- 


1.1984 


1.39 


1.44 


- 


0.03 


No 


383637 


Rl 


2.0 


1/30 


1.0074 


0.28 


1.25 


0.0018 


0.00 


Yes 


18021 


R2 


3.0 


1/30 


1.0224 


0.32 


1.17 


0.0025 


0.00 


Yes 


28449 


R3 


5.0 


1/30 


1.2241 


1.05 


1.21 


0.0197 


0.14 


No 


180004 


R4 


10 


1/30 


1.4773 


5.14 


1.72 


0.0578 


0.48 


No 


490026 


RTl 


2.0 


1/3 


1.0665 


0.42 


0.99 


0.0081 


0.01 


No 


174270 


RT2 


3.0 


1/3 


1.1537 


0.77 


1.22 


0.0080 


0.01 


No 


212730 



(Yoon & Langer "2002*1, overcritical rotation is even reached 
earlier. 

Livio & Pringle J1998> suggested nova explosions as a pos- 
sible mechanism for angular momentum loss to explain the 
white dwarf rotation velocity in typical CV systems. However, 
within the single degenerate SNe la progenitor scenario in 
which the white dwarf must grow to the Chandrasekhar limit 
by mass accretion, we need another mechanism to solve the 
angular momentum problem. To prevent the white dwarf from 
rotating over-critically, in our calculations, we limited the an- 
gular momentum gain in a simple way as follow. If the sur- 
face velocity at the white dwarf equator is below the Keplerian 
value, the specific angular momentum of the accreted matter 
is assumed to have the Keplerian value. Otherwise, no angu- 
lar momentum is allowed to be accreted onto the white dwarf 
while mass is still accreted: 



J (accreted matter) 



^Kcplcr^WD if Vs < WKcplc 
if Ws = WKcplo 



(1) 

where WKopicr = v^GA/wd/^wd, and Vg is the surface rota- 
tion velocity of the accreting white dwarf at its equator. This 
assumption may be justified by considerations of Paczyhski 
(1991 1 and Popham & Narayan ( 1991 1 to allow angular mo- 
mentum transport from the white dwarf surface through the 
boundary layer into the accretion disk by turbulent viscosity 
when the spin rate of the white dwarf reaches the Keplerian 
value, without preventing efficient mass accretion. 

As discussed in Heger et al. ("2000"), the contribution of the 
rotationally induced instabilities to the chemical diffusion is 
considered by introducing a factor /c G [0, 1], which describes 
the ratio of the chemical diffusion coefficient to turbulent vis- 
cosity. Heger et al. (2000,1 and Heger & Langer (2000 1 chose a 



value of /c = 1/30 following Chaboyer & Zahn (lT992t (See 
also Pinsonneault et al. il989> . In the present study, we consider 
two different values for /c, i.e., 1/30 and 1/3 (see Table [0, in 



order to investigate how the strength of chemical mixing affects 
the nuclear shell burning. 

For a comparative study, rotating and non-rotating models 
with the same initial conditions have been investigated. The 
properties of all computed models are summarized in Tabled 



4. Evolution 

The results of our calculations are summarized in Table[2 For 
sequences NRl, NR2, NR3, Rl and R2, we followed the evo- 
lution until the helium shell flashes become strong enough to 
drive the surface luminosity to the Eddington limit. For the 
rest, calculations are stopped after a large number of models 
has been calculated (cf. Tab. [3. 

4. 1. Helium shell burning and luminosity 

In Fig. n the nuclear luminosity due to the helium burning in 
the sequences with M = 2 x IQ-"^ Mg/yr (NRl, Rl & RTl) 
is shown as a function of the total mass of white dwarf. In 
sequence NRl, helium shell burning is not stable but under- 
goes thermal pulses, which are driven by the thin shell instabil- 
ity (e.g. Schwarzshild & Harm 1 19651 Weigert[T966j. It is wefl 
known that this phenomenon occurs also in intermediate mass 
stars in their thermally pulsating asymptotic giant branch (TP- 
AGB) phase (e.g. Iben & Renzini 1983 1. The thermal pulses 
in sequence NRl become stronger as the white dwarf accretes 
more matter (cf. Fujimoto & Sugimoto"1979"l. Finally, the shell 
flash at the sixth pulse drives the surface luminosity to the 
Eddington limit when AfwD — 1.002 M0 (cf. Fig.|2j, at which 
point we stop the calculation. Perhaps, the white dwarf would 
experience a radiation driven optically thick outflow at this 
stage, which would reduce the mass accumulation efficiency 
(Kato & Hachisu ri999> . It is beyond the purpose of this paper 
to consider the ensuing evolution in detail. 
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Fig. 1. Evolution of nuclear luminosity due to helium burn- 
ing in the helium accreting white dwarf models with M = 
2 X 10~^ Mo/yr. The box in panel (a) gives the result from 
the non-rotating case (NRl) while the results from the rotating 
models with f^ ~ 1/30 (Rl) and fc = 1/3 (RTl) are shown in 
panels (b) and (c). The abscissa denotes the total mass of white 
dwarf, which serves as a measure of time. 



In the rotating sequence Rl with the same accretion rate, 
thermal pulses become significantly weaker compared to the 
non-rotating case. For instance, the helium shell source reaches 
log Lj^/Lq ~ 6.5 in the non-rotating case when the white 
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Fig. 2. Evolution of white dwarf in the sequences NRl (a) and 
Rl (b) in the Hertzprung-Russell diagram. The starting point is 
marked by a filled circle. The ending point, where the surface 
luminosity reaches the Eddington limit, is marked by an aster- 
isk for each case. The ruggedness which appears in some parts 
of the evolutionary tracks is due to the limited number of digits 
in the numerical outputs. 



dwarf has accreted about 0.004 M©, while the rotating se- 
quence reaches only logL^/LQ ~ 4.0 at the same epoch 
(see Fig. Q. As a consequence, the white dwarf can accrete 
more mass (AMacc — 0.0096 Mq) than in sequence NRl 
(AMacc — 0.0039 M0) before the surface luminosity reaches 
the Eddington limit. In the HR diagram (Fig.|2}, the rotating se- 
quence is located in lower effective temperatures than the non- 
rotating sequence, since the white dwarf envelope is expanded 
by the centrifugal force. 

The helium shell source becomes more stable with stronger 
chemical mixing. Results of sequence RTl in Fig.[l];, where 
the chemical mixing is 10 times stronger than in sequence Rl, 
show that the helium is burned more or less steadily until Mwd 
~ 1.01 M0 . More than 1 00 pulses were followed thereafter, but 
the surface luminosity still did not reach the Eddington limit. 
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Fig. 3. Same as in Fig.[T]but with M = 3 x 10"'^ Mg/yr (NR2 
(a), R2 (b) and RT2 (c)). 



The evolution of sequence NR2 (M = 3 x 10"^ Mg/yr) 
shown in Fig.|5]is similar to that of NRl, only the white dwarf 
reaches the Eddington limit somewhat later, i.e., when AfwD — 
1 .004 Mq . The thermal pulses are weakened in the correspond- 
ing rotating sequence with /c — 1/30 (R2, Fig.lSji), and about 
4 times more mass could be accreted in the white dwarf of se- 
quence R2 before reaching the Eddington limit. With a stronger 
chemical mixing (fc — 1/3, RT2, Fig.|3j;), steady helium shell 
burning continues until AfwD — 1.01 M©, and weak thermal 
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Fig. 5. Evolution of white dwarf in sequence R4 in the 
Hertzsprung-Russell diagram. The region with logLs/i© = 
3.4 . . . 4.2 and logT = 5.46 . . . 5.48 represents the epoch from 
the onset of helium accretion to the first mild helium shell flash 
(cf. Fig. Hi). 



pulses with logLN/-^© ^ 5.0 follow. The evolution of this se- 
quence was stopped when Mwd — 1-15 Mq, as the number of 
thermal pulses before reaching the Chandrasekhar mass would 
exceed 10^. 

The stabilizing effect of rotation appears more prominently 
at M = 5 X 10^^ Mq/jt. In the non-rotating case (NR3), the 
helium shell burning is still unstable even with such a high ac- 
cretion rate, as shown in Fig.|^, and the white dwarf reaches 
the Eddington limit when Mwd — 1-013 M©. In the corre- 
sponding rotating sequence (R3), very weak thermal pulses 
appear when Mwd — 1-005 M©, but are suppressed when 
Mwd ^ 1-02 M©. Steady helium shell burning continues until 
A/wD — 1.14 M©, from which the shell source becomes again 
unstable. Note that the thermal pulses shown in Fig@J) look like 
noise but are actually many successive thermal pulses, each of 
which is resolved by about lO'^ models. 

With an accretion rate of M = 10~^ M©/yr, the non- 
rotating white dwarf undergoes unstable helium shell burning 
after A/wd ^ 1-03 M© (NR4, Fig.gJ;). We have foflowed the 
evolution until AfwD = 1-2 M© and find that thermal pulses 
do not become strong enough to lead the white dwarf to the 
Eddington limit. On the other hand, in the rotating sequence 
R4 (Fig.|4jl), the white dwarf can grow up to A/wd — 1.48Mq 
without encountering significant thermal pulses. 

The evolution of the white dwarf of sequence R4 in HR 
diagram is shown in Fig. |5l The surface luminosity remains 
nearly constant at logLg/L© ~ 4.2, while the surface temper- 
ature increases from logT — 5.46 to logT — 5.88, with which 
the white dwarf will appear as a super-soft X-ray source (e.g. 
Kahabka & van den Heuvel [T997l Greiner l2000t . 

Results presented in this section show that the helium shell 
source becomes more stable with rotation, and that the stronger 
the chemical mixing the more stable the shell burning becomes. 
We discuss the physical reasons for this behavior in Sect.|6| 
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Fig. 4. Evolution of nuclear luminosity due to helium burning in the helium accreting white dwarf models for A/ = 5 x 
10~'^ Mg/yr (NR3 (a) and R3 (b)) and M = lO"*^ Mg/yr (NR4 (c) and R4 (d)). The abscissa denotes the total mass of white 
dwarf, which serves as a measure of time. Note that what appears to look like noise in the plots for R3 and NR4 corresponds 
to many successive thermal pulses each of which being as well resolved in time (i.e., by roughly 10"^ time steps) as those in the 
plots in Fig.[l]and|3l 



4.2. White dwarf spin and angular momentum 
transport 



According to Yoon & Langer J2004a> . who considered ac- 
cretion of carbon-oxygen rich matter, non-magnetic accreting 
white dwarfs with M > 10^^ Mo/yr show the following 
features in relation to the angular momentum redistribution, 
a) Eddington-Sweet circulations, the secular shear instability, 
and the GSF instability are dominating the angular momen- 
tum redistribution in the non-degenerate outer envelope, while 
the dynamical shear instability is the most important in the de- 
generate core, b) The degree of differential rotation in the de- 
generate core remains strong throughout the evolution, with a 
shear strength near the threshold value for the onset of the dy- 
namical shear instability, c) The resulting differential rotation 
changes the white dwarf structure such that the center does not 
reach carbon ignition even when the white dwarf reaches the 
Chandrasekhar Hmit (-- 1.4 M©). d) More than 60 % of the an- 
gular momentum of the accreted matter is rejected by the con- 
dition that the white dwarf should not gain angular momentum 
when its surface rotates critically (Eq.[l). 



It is found that all these features remain the same in the 
present study, where helium accretion is considered. Fig. |6l 
shows angular velocity profiles throughout white dwarf mod- 
els of sequence RTl, R3 and R4, at various times, where the 
angular velocity is given in units of the local Keplerian value 
(i.e., cjKcp = \/GMr/r^). Note that the spin rate remains well 
below the local Keplerian value throughout the whole star in all 
cases. 

The white dwarf model with Mwd — 1.4 Mq in se- 
quence R4 has a central density of 3.54 xlO® g/cm'^, which 
is still far from the carbon ignition density, (cf. Ostriker & 
Bodenheimer [T968l Durisen [T975l Durisen & Imamura flMTt . 
As shown in the mentioned papers, a differentially rotating 
white dwarf can be dynamically stable up to ^ 4.0 Mq. As dis- 
cussed in Yoon & Langer J2004al l. this suggests that the secular 
instability to the gravitational wave radiation reaction, which 
is often called as the 'CFS instability' (Chandrasekhar 0*9701 



Friedman & Schutz [19781 1. might be important for the final fate 
of accreting white dwarfs at near or above the Chandrasekhar 
limit. The ratio of the rotational energy to the gravitational en- 
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Fig. 6. Spin velocity at different white dwarf masses, normal- 
ized to the local Keplerian value, as function of the mass coor- 
dinate for RTl, R3 and R4. 



ergy (i?rot/|£'grav|), which determines the growth time scale of 
the CFS instability via the bar-mode or the r-mode (see Yoon 
& Langer "2004a"), remains smaller than 0. 1 in the white dwarf 
models of sequence R4, and the last model (A/wd = 1-48 M0) 
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Fig. 7. Logarithm of the mass fractions of helium and carbon 
as function of the mass coordinate (solid line), in the white 
dwarf models of sequences NR4 (upper panel) and R4 (lower 
panel) at a time when AfwD = 1-0305 M©. The dotted and 
dashed lines denote the corresponding energy generation rate 
(right scale) for the 3a and ^^C(a, 7) reactions respectively. 



has £'rot/|£^grav| = 0.06. Sincc ii'rot / 1 -Egrav | ^ 0.1 is required 
in white dwarfs with strong differential rotation for the onset 
of the bar-mode instability (Imamura et al. 119951 . only the r- 
mode may be relevant at this point. The expected growth time 
of the r-mode instability in this model is about 1 . 1 x 10^ yr (see 
Yoon & Langer 2004a cf. Lindblom lT999t . We refer to Yoon 
& Langer (2004a i for a detailed discussion of the final fate of 
rapidly rotating massive white dwarfs. 



5. Chemical mixing and nucleosyntliesis 

The mixing of a-particles into the CO-core and of carbon into 
the helium envelope has consequences for the helium burning 
shell source. Here, we discuss the resulting changes in chem- 
ical structure of the white dwarf, while in the next section we 
analyze how rotation and rotationally induced chemical mixing 
affects the stability properties of the shell source. 
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Fig. 8. Carbon mass fraction as a function of the mass coor- 
dinate for the last computed models of sequences R4 (dashed 
line) and NR4 (dotted line). The thin solid line corresponds to 
the initial model. Note that the mass of the helium envelope is 
only of the order of 10"'* M0, which is not resolved in this fig- 
ure. The sharp peak in the carbon abundance near the surface 
at each white dwarf model is due to the active 3a reaction at 
the bottom of the helium envelope (cf. Fig.0. 



Secondly, the enhanced *^C(a,7)-efficiency results in an 
increased oxygen abundance in the ashes of the helium burn- 
ing shell, at the expense of carbon. Fig. |8l shows the carbon 
mass fraction throughout the white dwarf for the last computed 
models of sequences NR3 and R3. In the non-rotating model, 
the carbon mass fraction produced by the shell source reaches 
almost Xc — 0.68 at a white dwarf mass of 1.03 Mq. At 
this time, the helium shell source becomes unstable and devel- 
ops thermal pulses. As the layers at and above the shell source 
become convectively unstable during each pulse, the resulting 
carbon abundance is reduced. The convective mixing becomes 
more and more significant as thermal pulses become stronger 
(cf. Fig.|^), which is reflected in the drop of Xq from about 
0.68 to 0.55 in the range 1.03 Mq < M, < 1.20 Mq in Fig.Hl 
A comparison with the rotating model shows that the carbon 
mass fraction is mostly below 0.5, even though the rotating 
model does not suffer from thermal pulses for most of the time. 
While the average carbon mass fraction in the accreted layer 
is 0.58 in the non-rotating model, it is 0.45 in the layer with 
1.0 < A/^/Mq < 1.2 of the rotating model. The mean value 
of Xc throughout the accreted layer in the rotating model (i.e., 
1.0 < Mr/MQ < 1.48) is 0.38. 

The carbon and oxygen mass fractions in the accreted en- 
velope of a SN la progenitor are important from the point of 
view of nucleosynthesis and supernova spectroscopy, and may 
even affect the peak brightness of the supernova explosion (cf. 
Umeda et al. [T9^ Hoflich et al. l2000l Dominguez et al. l20Un . 



Fig-0shows the chemical structure and the nuclear energy 
generation rate in the helium burning shell in the white dwarf 
models of sequences R4 and NR4 when AfwD = 10305 Mq. 
The helium layer in the rotating model extends farther into the 
CO core, compared to the non-rotating case, due to the rotation- 
ally induced chemical mixing. This renders the ^^C(a, 7)^^0 
reaction more active in the rotating model, with two main con- 
sequences. 

Firstly, the contribution of the ^^C{a, 7)*^0 reaction to the 
total energy production is considerably increased. In the given 
models in Fig.0 the rotating model gives log Lsq/Lq = 3.938 
and log_Li2c(Q,-Y)/-LQ = 3.562, while the non-rotating model 
gives log Lsa/L© = 3.963 and log Li2c(a^T,)/i0 = 3.338, 
where L^a and Li2(j(^^.y-j are the integrated nuclear energy gen- 
eration rate over the shell source, due to the triple alpha reac- 
tion and the *^C(q;, 7)*^0 reaction, respectively. The 3a reac- 
tion is weaker in the rotating model than in the non-rotating 
one due to a lower density in the helium envelope (cf. Fig.|9jl 
and due to a smaller mass fraction of helium at the peak of the 
nuclear energy generation rate caused by rotationally induced 
mixing (Fig.Q. However, the increase of the energy genera- 
tion by the *^C(a,7)^^0 reaction is more important on the 
whole, leading to a higher nuclear luminosity in the rotating 
model (logi^N/LQ = 4.094) compared to the corresponding 
non-rotating case (logLN/L© = 4.060), for the models given 
in Fig.0 Furthermore, the shell source is hotter in the rotating 
model [Tshcii/(10^ K) ~ 2.16] than in the comparable non- 
rotating model [rshcii/(10^ K) = 2.03], which affects the sta- 
bility of the shell source as discussed in Sect.|6l 



6. Stability of hieiium slieii burning 

Our results show that helium shell burning is more stable with 
rotation, compared to the non-rotating case (Figs.[n|3]EJ. The 
stability of a helium shell source is mainly determined by three 
factors: geometrical thickness, degree of degeneracy and tem- 
perature of the shell source. In other words, a helium shell 
source is more stable if it is thicker, less degenerate and hot- 
ter. While a general and quantitative stability criterion for shell 
sources is worked out in a separate paper (Yoon, Langer & van 
der Sluys 2004"), we discuss the qualitative effects of rotation 
on the shell source stability of the models presented in Sect. 3 
in the following. 

Fig. m shows the evolution of the helium shell sources of 
sequences NR4 and R4 in the plane of the degeneracy param- 
eter T] (:= ip/kT, e.g. Clayton[l968) and the relative thickness 
of the shell source D/r^, where D and rg are the thickness 
and the radius of the shell source, respectively. Here, the en- 
ergy generation rate weighted mean over the shell source is 
used for the degeneracy parameter The thickness of the shell 
source is defined such that the energy generation rate at each 
boundary is 2 x 10"'^ times its peak value. Fig.|9]illustrates the 
evolution toward instability for the example of sequences NR4 
and R4. In sequence NR4, the mean temperature in the shell 
source remains nearly constant, around T ~ 2.02 x 10^ K dur- 
ing the stable shell burning phase. Fig.|9]shows the border be- 
tween stability and instability for T = 2 x 10^ K, as obtained 



by Yoon, Langer & van der Sluys J2004> . as dotted line. The 
shell source becomes more degenerate and thinner as the white 
dwarf mass grows, and finally enters the unstable regime when 
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Fig. 9. Evolution of the helium shell source in the sequences NR4 and R4 during stable shell burning phase, in the plane spanned 
by the degeneracy parameter 77 and the relative thickness of shell source D/rg- Here D denotes the thickness of shell source and 
Tg its radius. Each asterisk and black dot represents one selected stellar model in sequence NR4 and R4, respectively. The white 
dwarf mass in solar masses and the temperature of the shell source in the unit of 10^ K are indicated for various models. The thin 
lines separate the stable and unstable regime, for T/(10^ K) = 2.0 (dotted line), 2.2 (short dashed line), 2.4 (dotted-dashed line), 
2.6 (three dotted - dashed line) and 2.8 (long dashed line) according to the analytic stability criterion of Yoon, Langer & van der 
Sluys J2004I I. I.e., helium shell burning is predicted to be stable (unstable) above (below) the lines, while thermal pulses become 
significant from Af-wo — 1-03 M© in sequence NR4, and weak thermal pulses appear from Mwd — 1-41 M© in sequence R4. 



AfwD — 1.03 M0, resulting in the onset of thermal pulses as 
shown in Fig.|^. 

As already discussed in Sect. |5J the shell source in the 
rotating models extends over a wider mass range due to ro- 
tationally induced chemical mixing. I.e., significant amounts 
of a-particles are mixed into the CO core and participate in 
the ^^C(a,7)^^0 reaction. The presence of the centrifugal 
force lowers the density, and were it not for the effect of mix- 
ing, the shell source temperature would decrease with rotation. 
However, due to the enhanced ^^C(a, 7)^^0 reaction, temper- 
ature and nuclear energy production are increased, despite the 
lower density. In summary, rotation changes all three physi- 
cal properties which affect the stability — thickness, degener- 
acy and temperature — such that stability becomes more likely. 
Fig.|9]shows that this effect is quite significant. While the evo- 
lutionary tracks of sequences NR4 and R4 are close to each 
other in the shell thickness-degeneracy plane, at a given mass 
the shell thickness is much larger and the degeneracy much 
lower in the rotating models. For example, at AfwD = 1.01 Mq, 
the non-rotating model gives ry ~ —2.9 and D/ts — 0.225 



and the corresponding rotating model gives 77 ~ —3.0 and 
D/ts — 0.281. Furthermore, we have 77 ~ —2.8 and D/ts — 
0.19 at AfwD = 1-03 M0, from which thermal pulses become 
significant in the non-rotating sequence, while 77 ~ —2.9 and 
D /ts ~ 0.25 even when AfwD is as large as 1.14 M© in the 
rotating model (Fig.|9j. 

Additionally, the instability border which is given by the 
thin lines in Fig.|9lis located lower at a given mass in the ro- 
tating case, due to the higher shell source temperature in the 
rotating models, which leads to a later onset of the instability 
even for a given shell source thickness and degeneracy. For in- 
stance, the shell source temperature in the rotating model when 
Mwd = 1-01 Mq is about 2.1 x 10^ K, which is higher than in 
the non-rotating models, where Tghcii remains at ~ 2.0 x 10^ K 
during the steady shell burning phase. The last analyzed model 
in sequence NR4 is already in the unstable regime (i.e., below 
the dotted line; D /vg — 0.19) and thermal pulses become sig- 
nificant right after as expected from the stability criterion. On 
the other hand, in sequence R4, the shell source temperature 
increases up to 2.6 x 10* K when A/wd = 1.41 M©, which 
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explains why the shell source is stable until ^fwD — 1-41 M0 
even though D /vg at this moment is smaller than in the last an- 
alyzed model of NR4. As shown in Fig.|^ (see also Fig.|5), the 
ensuing thermal pulses which appear when Mwd ^ 1-41 Mq 
in sequence R4 remain very weak until the last model (Mwd 
= 1.48 M© with Tshcii ^ 2.6 x lO'^ K), which is still near the 
stabiHty border line for T = 2.6 x 10^ K in Fig.|9l 

In summary, rotation stabilizes the instability of the helium 
shell burning in our models significantly. The shell sources in 
realistic white dwarfs, however, may be even more stable than 
our models predict, due to two reasons. Firstly, as outlined in 
Sect. 12] our method forces us to limit the structural effects of 
rotation to that obtained at about 60% of critical rotation. This 
means that we underestimate the centrifugal force in the layers 
which contain the shell source. Secondly, our models are one- 
dimensional, which restricts the pressure response due to any 
perturbation to the radial direction — which is generally not 
true in the realistic situation. In a rapidly rotating star, the pres- 
sure can also respond to a perturbation horizontally, unless the 
shell source is located at the poles. We believe that this multi- 
dimensional effects can stabilize the shell source further. 

7. Final discussion 

We have performed numerical simulations of helium accret- 
ing CO white dwarfs with an initial mass of 0.998 M© and 4 
different mass accretion rates (2, 3, 5 and 10xlO~^ MQ/yr), 
considering spin-up by mass accretion. The main feature of the 
spin-up is the angular momentum transport from the accreted 
matter into the white dwarf, which creates differential rotation 
in the white dwarf interior. This induces transport of chemical 
species from the accreted matter into the CO core, through the 
shear instability, Eddington-Sweet circulations, and the GSF 
instability. 

It is shown that the helium burning layer becomes spatially 
more extended compared to the non-rotating case, since the ac- 
creted a-particles are mixed into the carbon enriched layer, par- 
ticipating actively in the ^^C(a, 7)^^0 reaction (Sect.|5|l. The 
carbon-to-oxygen ratio is significantly reduced in the accreted 
layer The centrifugal force lowers density in the shell source 
and degeneracy is significantly lifted. The total energy gener- 
ation in the helium shell source increases with rotation, due 
to the increased contribution from the ^^C(q;, 7)^^0 reaction, 
rendering the helium shell source as hot as in the non-rotating 
cases, despite the decreased density. All these effects contribute 
to stabilizing the helium shell burning significantly (Sect. |6j. 
Thermal pulses, if they occur, are considerably weakened with 
rotation. 

As mentioned previously, and as discussed by Cassisi et 



al. ([1998), Kato & Hachisu ( fT999t and Langer et al. ( I2002t . a 
high mass accumulation efficiency on a white dwarf by hydro- 
gen accretion has been seriously questioned due to the unstable 
helium shell burning which is found when helium is accreted 
at a rate which gives steady hydrogen shell burning. Note that 
the accretion rate of 5 x 10^^ M0/yr corresponds roughly to 
the maximum rate for the steady hydrogen burning in a 1 .0 M© 
white dwarf (Nomoto'1982a"). Our results indicate that even if 
hydrogen is accreted and burned steadily, the subsequent he- 



lium shell burning may be unstable, with strong helium shell 
flashes driving the stellar envelope to the Eddington limit, un- 
less effects of rotation are considered. Kato & Hachisu ( I1999> 
suggested that if the new OPAL opacities are adopted, an opti- 
cally thick wind can remove mass from the system without the 
white dwarf filling its Roche lobe. They also suggested that the 
effective mass accumulation rate can remain high enough for 
the white dwarf to grow to the Chandrasekhar limit. However, 
their results are based on a very massive white dwarf model of 
L3 M0. Cassisi et al. (' 1998t . on the other hand, showed that 
hydrogen accreting white dwarfs with initial masses of 0.56 
M0 and 0.8 Mq fill their Roche lobe once helium flashes occur, 
leaving doubts about the efficient mass increase of the white 
dwarf from a relatively low mass by hydrogen accretion. 

Our results imply that rotation may help white dwarfs to 
grow in mass by stabilizing helium shell burning. In the canon- 
ical single degenerate scenario, however, a main sequence star 
or a red giant is supposed to transfer hydrogen rich matter 
onto the white dwarf (Li & van den Heuvel [19971 Hachisu 



et al. 1999; Langer et al. I2000t . We still need to investigate 
the case of hydrogen accretion, in order to answer the ques- 
tion whether rotation may be a potential solution in explaining 
the observed SNe la rate from the considered binary systems. 
Since the stability conditions for hydrogen shell burning are not 
qualitatively different from that for helium shell burning (Yoon, 
Langer & van der Sluys 2004 1, we expect that the effects of ro- 
tation on the hydrogen shell source are likely to give a similar 
conclusion as for the helium shell burning. 

On the other hand, hydrogen accretion with a rate of 10"^ 
Mo/yr is supposed to expand a 1.0 Mq white dwarf into a red 
giant phase (Nomoto ll982al Nomoto & Kondo 1991, Cassisi 
et al. 1998). In a close binary system, the hydrogen envelope 
expanded to the red giant phase will fill the Roche lobe, caus- 
ing significant mass loss or the merging of both stars, unless the 
mass accretion rate is not reduced by an optically thick stellar 
wind, as suggested by Hachisu et al (1996', I1999t . However, 
steady helium accretion with such a high accretion rate can be 
realized, for example, in a binary system which consists of a 
white dwarf and a helium giant star ("helium Algols", Iben & 
Tutukov ,1994j . in which the mass transfer rate can amount to 
IQ-^ - IQ-^ Me/yr, as shown by Yoon & Langer ([20031 . Our 
results indicate clearly that the mass accumulation efficiency 
can be significantly enhanced by rotation for helium giant star 
-I- white dwarf binary systems. Although the calculated produc- 
tion rate of SNe la via this route is currently rather uncertain 
(Iben & Tutukov [T994[ Branch et al. 1998) . consideration of 
rotation may increase it strongly. 

In the present study, we did not consider stellar winds 
from white dwarfs, on which rotation might cause significant 
effects. When the stellar surface approaches the critical ve- 
locity, the mass loss due to radiation-driven winds could in- 



crease dramatically (Friend & Abbott 1986 Langer I1997> . 
Since stellar winds may also carry a large amount of angu- 
lar momentum (Langer "1998), this would involve a compli- 
cated history of angular momentum gain or loss in accreting 
white dwarfs: accretion will spin up white dwarfs, while stel- 
lar winds will brake them. It is also possible that the physi- 
cal conditions in the shell source in a rapidly rotating white 
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dwarf vary with latitude, which may lead to interesting con- 
sequences such as non-spherical stellar winds (cf. Owocki & 
Gayely 1997 Maeder 1999 1. We will investigate these effects 
systematically in near future. 

Finally, we note that an increase in the uncertain efficiency 
of rotationally induced chemical mixing results in more sta- 
ble shell burning (Fig. [2 and |3}, than the case where the cal- 
ibration by Heger et al. ©00 1 and Heger & Langer fiOOOit 
is adopted. Their calibration gives a similar chemical mixing 
efficiency as in Pinsonneault et al. ( 1989 1 who made a calibra- 
tion to reproduce the solar surface ^Li abundance, as well as 
to the theoretical estimate by Chaboyer & Zahn (1992i. The 
physical conditions in the helium shell source in white dwarfs 
do not qualitatively differ from those in main sequence stars. 
However, for accreting white dwarfs, another mixing mecha- 
nism which was not considered in the present study may be im- 
portant: magnetic instability induced by differential rotation in 
radiative layers (Spruit 2002 1. Heger et al. J2003t and Maeder 



& Meynet (I2003> show that the turbulent diffusion can become 
stronger by several orders of magnitude due to the magnetic 
instability, compared to the case where only rotational effects 
are considered. The differential rotation induced magnetic in- 
stability can be even more important in accreting white dwarfs, 
in which a strong shear motion is retained due to the continu- 
ous angular momentum transport from the accreted matter, as 
shown in Sect. 14.21 The efficiency of the chemical mixing may 
increase accordingly, which may stabilize the shell source even 
further than in the present study. This possibility is currently 
under investigation. 
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